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We present results from a weak lensing analysis of the galaxy cluster A1689 (z = 0.183) 
based on deep wide-field imaging data taken with Suprime-Cam on Subaru telescope. A 
maximum entropy method has been used to reconstruct directly the projected mass dis- 
tribution of A1689 from combined lensing distortion and magnification measurements of 
red background galaxies. The resulting mass distribution is clearly concentrated around 
the cD galaxy, and mass and light in the cluster are similarly distributed in terms of shape 
and orientation. The azimuthally-averaged mass profile from the two-dimensional recon- 
struction is in good agreement with the earlier results from the Subaru one-dimensional 
analysis of the weak lensing data, supporting the assumption of quasi-circular symmetry 
in the projected mass distribution of the cluster. 
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1. Introduction 

Weak gravitational lensing of background galaxies provides a unique, direct way 
to study the mass distribution of galaxy clusters. Recent improvements in the 
quality of observational data usable for lensing studies now allow an accurate de- 
termination of themass distribution in clusters. A1689 is one of the best studied 
lensing clusters j^^^^EE] located at a moderately low redshift oi z = 0.183. Deep 
HST/ACS imaging of the central region of A1689 has revealed 106 multiply lensed 
images of 30 background galaxies, which allowed a detailed reconstruction of the 
mass distribution in the cluster core (10/i^^kpc ^ 200ft-^^kpc).^ In Ref. jH we 
developed a model-independent method for reconstructing the cluster mass pro- 
file using azimuthally-averaged weak-lensing distortion and magnification measure- 
ments, and derived a projected mass profile of A1689 out to the cluster virial radius 
[r ;S 2ft,~^ Mpc) based on the wide-field, deep imaging data taken with Suprime- 
Cam on the 8.2m Subaru telescope. The combined strong and weak lensing mass 
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profile is well fitted by an profile with high concentration of Cvir ~ 13.7, 

which is significantly larger than theoretically expected (Cvir — 4) for the standard 
LCDM modelP 

In this paper we present a weak lensing analysis of A1689 using wide-field Sub- 
aru imaging data, with special attention to the map-making process. Throughout 
this paper, we use the AB magnitude system, and adopt the concordance ACDM 
cosmology with (Jlnio = 0.3, fi^o — 0.7, h — 0.7). In this cosmology one arcminute 
corresponds to the physical scale 129kpc//i for this cluster. 

2. Sample selection 

For our weak lensing analysis we used Subaru/Suprimc-Cam imaging data of A1689 
in V (1,920s) and SDSS i' (2,640s) retrieved from the Subaru archive, SMOKA 
(see Ref. HEIfor more details). The FWHM in the final co-added image is 0"82 
in V and 0"88 in i' with 0"202 pix-\ covering a field of 30' x 25'. The Hmiting 
magnitudes are V = 26.5 and i' = 25.9 for a Scr detection wit hin a 2"aperture. A 
careful background selection is critical for a weak lensing analysis For the number 
counts to measure magnification, we define a sample of 8,907 galaxies (fi^ = 12.0 
arcmin^^) with V — i' > 1.0. For distortion measurement, we define a sample of 
5,729 galaxies {fig = 7.59 arcmin"^) with colors 0.22 mag redder than the color- 
magnitude sequence of cluster E/SO galaxies, {V -i') + 0.0209i' -1.255 > 0.22. The 
smaller sample is due to the fact that distortion analysis requires galaxies used are 
well resolved to make reliable shape measurement. We adopt a limit of i' < 25.5 to 
avoid incompleteness effect . In what follows we will assume (zg) = 1 for the mean 
redshift of the red galaxies ,^3131 but note that the low redshift of A1689 means that 
for lensing work, a precise knowledge of this redshift is not critical. 

3. Lensing Distortions 

We use the IMCAT package developed by N. Kaiser @ to perform object detection, 
photometry and shape measurements, following the formalism outlined in Ref.i9j We 
have modified the method somewhat following the procedures described in Ref. [TOl 
To obtain an estimate of the reduced shear, ga = 7q/ (1 — k), we measure the image 
ellipticity from the weighted quadrupole moments of the surface brightness of 
individual galaxies. Firstly the PSF anisotropy needs to be corrected using the star 
images as references: 

e'^^e^- PZ^q} (1) 

where Psm is the smear polarizability tensor being close to diagonal, and g* = 
{P*sm)~^^e.i is the stellar anisotropy kernel. We select bright, unsaturated foreground 
stars identified in a branch of the half-light radius {rh) vs. magnitude ii') diagram 
(20 <i' < 22.5, (r-,0„,edian 2.38 pixels) to calculate g^. 

^http://www.ifa.hawaii/kaiser/IMCAT 
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Fig. 1. Gravitational reduced-shear field in A1689 obtained from shape distortions of the red 
background galaxies, smoothed with a Gaussian with FWHM = 2' for visualization purposes. The 
shaded circle indicates the FWHM of the Gaussian. The coordinate origin is at the optical center 
(see Figure |4]l. 



In order to obtain a smooth map of g* which is used in equation ([T]) , we divided 
the 9K X 7.4K image into 5x4 chunks each with 1.8K x 1.85K pixels, and then 
fitted the q* in each chunk independently with second-order bi-polynomials, q'^{0), 
in conjunction with iterative cr-clipping rejection on each component of the residual 
e* — P*smqfj{d)- The final stellar sample consists of 540 stars, or the mean surface 
number density of n» = 0.72 arcmin^^. From the rest of the object catalog, we 
select objects with 2.4 < r;i < 15 pixels as an i'-selected weak lensing galaxy 
sample, which contains 61,115 galaxies or fig ~ 81 arcmin"^. It is worth noting 
that the mean stellar ellipticity before correction is (ei*,e2*) ~ (—0.013,-0.018) 
over the data field, while the residual e* after correction is reduced to e"'''' = 
(0.47 ± 1.32) X 10-'', = (0.54 ± 0.94) x 10"^. The mean offset from the null 
expectation is |e*"^| = (0.71 ± 1.12) x 10^''. On the other hand, the rms value of 
stellar ellipticities, Ce* = (|e*P), is reduced from 2.64% to 0.38% when applying 
the anisotropic PSF correction. Second, we need to correct the isotropic smearing 
effect on image ellipticities caused by seeing and the window function used for the 
shape measurements. The pre-seeing reduced shear can be estimated from 

9c. - {Pg-')o.0e'^ (2) 

with the pre-seeing shear polarizability tensor P^^. We follow the procedure de- 
scribed in Ref. JJlto measure P^. We adopt the scalar correction scheme 
namely, P^^ = ^tT[P^]Sai3 = P^Sap- The Pg measured for individual objects are 
still noisy especially for small and faint objects. We remove from the galaxy cat- 
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alog those objects that yield a negative value of P| estimate to avoid noisy shear 
estimates. We then adopt a smoothing scheme in object parameter spacet^ ^*^ ! ^^ ! 
We first identify thirty neighbors for each object in Vg-i' parameter space. We then 
calculate over the local ensemble the median value 



Pg) of Pg and the variance 



o-g of g 
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ig2 using equation ([2]). The dispersion Ug is used as an rms error 
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of the shear estimate for individual galaxies. The mean variance Gg over the red 



galaxy sample is obtained as ~ 0.133, or ~ 0.36. Finally, we use the following 

estimator for the reduced shear: ga = e^/ {Pf}- 

For map-making, we then pixelize the distortion data into a regular grid of 
-^pix = 21 X 17 independent pixels, covering a field of w 30' x 24'. The pixel size 
is Apix — 1.'4, and the mean galaxy counts per pixel is ~ 15. The bin-averaged 
reduced shear is given as ga,i = ga{di) = SfeeceiH "fe^a.fe/ Z^fesceiH '^fe (a = l,2;i = 
1, 2, A^pix), where ga.fe is the estimate of the ath component of the reduced shear 
for the kth galaxy, and Uk — 1/(0"^ fc+a^) is its inverse- variance weight softened with 
a constant a. Here we choose a = 0.4.^131 In Figure[T]we show the reduced-shear field 
obtained from the red galaxy sample, where for visualization purposes the ga{Q) is 
resampled on to a finer grid and smoothed with a Gaussian with FWHM — 2'. 
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Fig. 2. Distribution of the lensing magnification bias n/no measured from a red-galaxy sample in 
the background of A1689, smoothed with a Gaussian with FWHM = 2' for visualization purposes. 
The shaded circle indicates the FWHM of the Gaussian. 



4. Magnification Bias 

Lensing magnification, ^Ji{6), influences the observed surface density of background 
galaxies, expanding the area of sky, and enhancing the flux of galaxies. In the sub- 
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critical regime,'^ "^^ the magnification fj, is given by = 1/ [(1 — k)^ ^ li ~ tI] • The 
count-in-cell statistics are measured from the flux-limited red galaxy sample (see 
Sj2l) on the same grid as the distortion data: Ni = N{< nicut', di) = J2kecciu ^ ^'^^^ 
mcut being the magnitude cutoff corresponding to the flux- limit. The normalization 
and slope of the unlensed number counts Nq(< mcut) for our red galaxy sample 
are reliably estimated as n^^ = 12.6 ± 0.23 arcmin^^ and s = dlog NQ{m) / dm — 
0.22 ±0.03 from the outer region > lO'.^lThe slope is less than the lensing invariant 
slope, s = 0.4, and hence a net deficit of background galaxies is expected: Ni/No = 
n^^i/rififi = /if^""^. The masking effect due to bright cluster galaxies is properly 
taken into account and corrected for.^ Figure [5] shows a clear depletion of the red 
galaxy counts in the central, high-density region of the cluster. Note we have ignored 
the intrinsic clustering of background galaxies, which seems a good approximation,^! 
though some variance is apparent in the spatial distribution of red galaxies. 



5. Two-Dimensional Mass Reconstruction 

The relation between distortion and convergence is non-local, and the co nver gence 
K derived from distortion data alone suffers from a mass sheet degeneracy.^^^ How- 
ever, by combining the distortion and magnification measurements the convergence 
can be obtained unambiguously with the correct mass normalization. Here we com- 
bine pixelized distortion and magnification data of the red background galaxies, 
and reconstruct the two-dimensional (2D) distribution of k using a maximum en- 
tropy method (ME M) ext ended to account for positive/negative distributions of 
the underlying field.^l^U3 ■\Ye take into account the non-linear, but sub-critical, 
regime of the lensing properties, k, and 7^. We take ~ n{di) as the image 
to be reconstructed, and express a set of discretized k- values as p = {k^} 
[i = 1, 2, A^pix)- The total log-likehhood function, F(p) = — ln>C(p), is expressed 
as a linear sum of the shear/magnification data log-likelihood&l^ and the entropy 
ternPl 

F(p) = /g(p)-HZ^(p)-a5(p,m), (3) 
/.-"ln/:.-EE^^^^^^I^' (4) 

2—1 

where .ga,j(p) and A'i(p) are the theoretical expectations for ga,i and Ni, respec- 
tively, Ug^i = <7g{9i) is the rms error for gi = gi{9i) + ig2{0i), and S'(p,m) is the 
cross entropy function for the positive/negative fieldp^the m is a set of the model 
parameters and a{> 0) is the regularization constant. The maximum likelihood so- 
lution, p, is obtained by minimizing the function i^(p) with respect to p for given 
a and m. We take rrii = const = m, and determine by iteration the Bayesian value 
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of a for a given value of m. We found that the maximum-hkehhood solution for the 
Bayesian a is insensitive to the choic e of m . In the following we set m to be 0.2. 
We note that the adopted MEM prioA^'^ ensures Hi —>■ in the noise-dominated 
regime, i^(p) ^ — aS'(p, m) (i.e., maximizing the entropy alone). In order to quan- 
tify the errors on the mass reconstruction we evaluate the Hessian matrix of the 
function F{p) at p = p, iJy (p) = §^r3^lp=P' fro™ which the covariance matrix 
of the parameters p is given by Cij = {AniAnj) = (i/)~^(p). Figure [3] displays 
the K map reconstructed with the MEM method. In Figure [4] we compare the con- 
tours of the reconstructed k (thick) and the z'-band luminosity density of the cluster 
sequence galaxies (thin) superposed on the z'-band image of the central region of 
A1689. 




X (arcmin) 

Fig. 3. The projected mass distribution k of A1689 reconstructed by a maximum entropy method 
using weak-Iensing distortion and magnification data of red background galaxies. For visualization 
purposes the K-map is resampled on to a finer grid, and smoothed with a Gaussian with FWHM = 
2'. The lowest contour is k = 0.07, and the contour steps are Are = 0.05. 



6. Model-Independent Mass Profile of A1689 

We show in Figure [5] the azimuthally-averaged profile of the reconstructed conver- 
gence, k{9), as a function of projected radius 9 from the optical center of A1689 
(circles). The vertical error bars represent the la uncertainties based on the error 
covariance matrix Cij of the reconstruction. Note that the error bars are correlated. 
Also shown for comparison are the earlier results from the HST/ACS strong lensing 
analysis (triangles) and the Subaru weak lensing analysis (squares) with the one- 
dimensional (ID) reconstruction method, respectively, along with the best-fitting 
model (solid) for the combined ACS-fSubaru profile (see Ref. [3|). 
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Fig. 4. Contours (thick) of the reconstructed lensing convergence k superposed on the i'-band 
image of the central region of A1689. The field size is 15' on a side. The lowest contour and the 
contour intervals are 0.05. Also superposed are the contours (thin) of the projected luminosity 
density of the cluster member galaxies defined by the color-magnitude sequence. North is to the 
top and East is to the left. 



7. Discussion and Conclusions 

We presented results from our weak lensing analysis of A1689 based on deep wide- 
field imaging data taken with Subaru/Suprime-Cam. We used a MEM algorithm to 
reconstruct the projected mass map in A1689 from combined distortion and mag- 
nification data of our red background galaxy sample. The combination of distortion 
and magnification data breaks the mass sheet degeneracy inherent in all recon- 
struction methods based on distortion information alone. Our results show that 
mass and light in A1689 are similarly distributed in terms of shape and orientation, 
and clearly concentrated around the cD galaxy (see Figure [4|). The resulting mass 
profile from the present full 2D reconstruction is in good agreement with the results 
from the earlier Subaru ID analysiJ^ (see Figure [5]), supporting the assumption of 
quasi-circular symmetry in the projected mass distribution. 
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Fig. 5. Reconstructed mass profile of A1689. The filled circles with error bars represent the results 
based on the 2D k map reconstructed with the MEM algorithm. The error bars are correlated. 
The triangle and square symbols with error bars show the results from the HST/ACS strong 
lensing analysis (Broadhurst et al. 2005, ApJ, 621, 53) and the Subaru weak lensing analysis 
with the ID reconstruction method (Broadhurst, Takada, Umetsu et al. 2005, ApJ, 619, L143, 
hereafter BTU05) respectively. The solid curve shows the best-fitting NFW profile for the combined 
ACS+Subaru data by BTU05 (M^ir = 1.93 x lO^^Mo.Cvir = 13.7). The best-fitting NFW profile 
for the ACS-|-Subaru profile has a high concentration, Cvir = 13.7, and somewhat overestimates 
the inner slope and is a bit shallower than the ID-based results at large radius. 
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